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We propose a versatile strategy for numerical renormalization group solution of general channel- 
mixing Kondo and Anderson models beyond previous reach, opening the door toward broad applica¬ 
tions in protocol non-perturbative machineries, such as dynamical cluster approximation and cluster 
dynamical mean field theory, for strongly correlated electron systems. We illustrate the strategy by 
investigating the quantum phase transitions in two quantum impurity models with cases untouched 
before. 
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The Kondo impurity model and its associated Ander¬ 
son impurity model are challenging many-body prob¬ 
lems and play an important role in condensed matter 
physics. [1] Such models develop infrared logarithmic di¬ 
vergence in perturbation theory, and hence require non- 
perturbative solutions. The Kondo impurity model was 
first solved by Wilson using the simultaneously invented 
numerical renormalization group (NRG). [2] The NRG 
captures the physics at exponentially decreasing energy 
scales iteratively, a key ingredient behind its great suc¬ 
cess. It proves to be one of the most accurate and ef¬ 
ficient methods for the quantum impurity problems, [3] 
and therefore has been broadly used as an impurity solver 
in the dynamical mean field theory (DMFT) and dynam¬ 
ical cluster approximation (DCA). [4-7] 

The quantum impurity is coupled to a noninteracting 
conduction band (or bath), the effect of which is to cause 
a so-called hybridization function V{uj). In the case of a 
scalar or diagonal matrix function T>(u>), a standard NRG 
procedure has been developed. [8-11] The key ingredient 
is to map the effect of T>(ui) into an open Wilson chain. 
However, the extension to quantum impurity models with 
■D(w) non-diagonal in the spin and/or local orbital basis is 
not straightforward at all. Such a matrix function is said 
to be channel-mixing henceforth, with a channel referring 
to a combination of spin and orbital (as well as atomic 
site in a cluster). This situation naturally arises in the 
presence of spin-orbital coupling, Cooper pairing, as well 
as in a cluster impurity. The intertwining spin, charge 
and orbital degrees of freedom are promising for harvest¬ 
ing novel quantum effects, and the cluster impurity is 
invoked in DCA and cluster DMFT (cDMFT). These in¬ 
teresting and important cases are barely addressed by 
NRG in the literature so far because of the lack of a ver¬ 
satile scheme to map a channel-mixing T>(u>) to a Wilson 
chain. The only exception to our knowledge, in fact a 
special one of the general cases, is when T>(ui) can be di¬ 
agonalized by a frequency-independent unitary transfor¬ 
mation in the presence of particular symmetries. [12-19] 
We also notice that the channel-mixing model we will be 
addressing is different to the usual multi-channel model 
where 2?(w) is actually channel-diagonal, although the 


latter is interesting on its own right. [20] 

In this Letter, we fill the gap caused by the above dif¬ 
ficulties. We propose a versatile scheme to map a gen¬ 
eral (matrix) hybridization function into an open Wil¬ 
son chain. We benchmark the mapping scheme against 
a nontrivial previous result, [12, 13] and we illustrate the 
scheme for two channel-mixing cases untouched before. 
Our strategy enables NRG solution of general channel¬ 
mixing Anderson and Kondo models, as would be desir¬ 
able in important protocols for strongly correlated elec¬ 
trons, such as DCA and cDMFT. 

A versatile mapping scheme : For definiteness we con¬ 
sider a generalized Anderson impurity model, while the 
Kondo impurity model will be addressed in the closing 
section. The Hamiltonian H = Hi mp + Ht, + Hhyb is 
composed of, 

Him P = Y,fXc fh + Hint, 
a,(3 

H b =J2 c laK b CKb, ( 1 ) 

k ,a,b 

Hhyb = Vn^ ^« 7 ^ aCk ’ a + hx -)- 

v k,a:,a 

Here f a is an annihilation field at channel a of the im¬ 
purity, Ck, a an annihilation field at channel a and at mo¬ 
mentum k in the bath, and N is the volume (or number 
of unit cells) in the bath. (The channel numbers in the 
impurity and the bath can differ in general.) As we men¬ 
tioned, we take a channel index as a combined label of 
spin and orbital (as well as the site in a cluster impu¬ 
rity). The concrete expressions for the matrices hi OCl hk 
and 7k, as well as the interaction Hi nt on the impurity 
are unnecessary at this stage. 

Integrating out the free bath leads to a (matrix) self¬ 
energy correction to the impurity, 

s (a) = X 7kG k (2)7k< ( 2 ) 

1 k 

where z = w ± i0 + in the retarded/advanced case, and 
Gk(z) = l/(z—hk) is the Green’s function at the complex 
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frequency z. Via Kramers-Kronig relation, S(z) can be 
completely characterized by 

V(uj) = i[ E(w + 10+) - S(w - i0+)]/27r. (3) 

which we call a (matrix) hybridization function. (For 
later convenience our definition differs to the usual one by 
a factor of 7 r.) However, 2?(w) is in general non-diagonal, 
i.e., channel mixing. Our purpose is to map the effect of 
such a 2?(w) into that of a generalized open Wilson chain. 

By definition, T>(uj) is hermitian and positive semi- 
definite. This enables us to write 

= l n > w )Pn(w)(n,w|, (4) 

n 

where p n (uj ) > 0 and \n,co) are the n-th eigenvalue and 
eigenvector of 2?(w). Similarly to the usual one-band 
case, [9, 10] we can re-express p n (u) as 



where e ra (x) is a continuous function of x for each n, 
subject to the requirement on the integration measures 

t 2 n (x) \dx\ = p n [e n (x)\ \de n (x)\. ( 6 ) 

[The modulus symbol is necessary if de n (x)/dx < 0.] The 
parametrization of e n {x) and t n (x) for p n (uj) is similar to 
the usual logarithmic way, [3, 9, 10] and is also provided 
in the Supplementary Materials. 

Substituting the expression of p n (uj) into V{w), we 
write 

T>(u)= dx y \n,u)tn(x)6[uj - e n {x)\{n,uj\ 

^ n 

= dx y |n, x)tn(x)6[u} - e n (x)\(n, x|, 

n 

where we made a replacement |n,w) —► | n,x) = 

|n, e n (x)}, valid in the presence of the delta-function in 
the integrand. This is a crucial step for the following 
discussions. In components, we have 

= / dx y V an (x)6[oj - e n (x)]V ra t / 3 (x), 

n 

where V an (x) = t n (x)(a\n,x). The above parametriza¬ 
tion maps Hb and Hh y b to, respectively, 

H b= dxyipl{x)e n (x)ip n (x ), 

J n 

H hyb= / dxy flv an {x)ijj n {x) + h.c., 

a,n 

where ip's are fermion fields in x. The mapping is ex¬ 
act since integrating out ip's leads to the same E(z) we 
started from. In practice, however, one would like to pro¬ 
ceed in a discrete space of x. Fortunately, as in the usual 
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FIG. 1. Graphical illustrations of (a) the discretized bath 
coupled to the impurity, (b) a channel-diagonal Wilson chain 
in the usual case, and (c) a channel-mixing Wilson chain in 
general cases. The number of channels is set as (but not 
limited to) 1 = 2 . 


case, the above mapping is designed such that its dis¬ 
cretized version provides a good approximation. For ex¬ 
ample, one can take Xj — (\j \ +ry)sign(j)[9] with nonzero 
integer j and a real twisting factor rj £ ( 0 , 1 ], and re¬ 
place the integration over x to a summation over j. Each 
case of 77 represents an approximation of the continuum 
model, and as far as 2 ?(w) is concerned, the average over p 
reproduces exactly the result from the continuum model. 
A graphical representation of a discrete mapping is shown 
in Fig. 1(a). 

We proceed to map the discretized model to an open 
Wilson chain, a key ingredient of NRG to reduce the com¬ 
putational cost and optimize the scaling behavior. In the 
special case when |n,x) = |n) is independent of x, a uni¬ 
tary transformation J2 a fa( a \ n ) fn makes H' b + H' hyb 
diagonal in n (while Hi mp may not and actually does 
not have to be so), a case accessible by the conventional 
mapping scheme. [12-19] But this does not apply to more 
general occasions we are aiming at. Let us assume the 
dimension of the single-particle Hilbert space of the im¬ 
purity is I (the number of channels), and that for H b is 
B = 2J/, where J is the number of Xj’s for j > 0 (or 
j < 0) retained in the discritization. Let us rewrite H b 
and H’ hyb compactly as 

H' hyb = F^ + ^V^F, 
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where 'I' and F are spinors composed of the if- and f- 
fields, respectively. In this form, £bxb and V/ X s are 
matrices. We perform a QR decomposition lA = UqTq, 
where TjJ is an upper triangular matrix. We use the 
columns of Uq as the first set of Krylov basis vectors 
to transform £, by block-Lanczos,[21, 22] into a block- 
tridiagonal matrix Et (with block size I) such that 
£ = U'EtU, where the columns of U' are Krylov ba¬ 
sis vectors, and the leading ones are from Uq. In the 
resulting Krylov space, or upon a canonical transforma¬ 
tion U'f! —»• $, we have 


K 


H' b + H' hyh ->Y, [$lE k <f> k + + h.c.], (7) 

fe=l 


where K is the number of block-Lanczos iterations, 
Q k denotes an /-component spinor such that $1 = 
(<I>{, <f> 2 , ‘ ‘ ‘ and we set $>o = F for brevity. For 

k > 1, E k (Xfc) is the /c-th / x I block element of Et 
along the diagonal (upper sub-diagonal). Notice that all 
Tfc’s are triangular matrices themselves. Eq. (7) defines 
the open Wilson chain we were after. For comparison, 
the usual Wilson chain is illustrated in Fig. 1(b), which is 
channel diagonal, while the Wilson chain in our general 
cases is illustrated in (c). In practice, the block-Lanczos 
procedure requires infinitely high precision and is trun¬ 
cated at a suitable stage K ~ J. Numerical examples of 
our mapping scheme for models discussed below can be 
found in the Supplementary Materials. 

Quantum impurity in an s-wave superconducting bath: 
Given the above versatile mapping, we are able to study 
any quantum impurity models, with or without channel 
mixing. Here we consider an Anderson impurity coupled 
to a conventional s-wave superconductor. The Hamilto¬ 
nian in the Nambu space is composed of 


H imp = lf'_f€f(T z 1pf - ^U{if^a z ipf) 2 , 
Hb = ^ V’k( £ kO-~ + AkCrJ'i/’k, 

k 


H hyb = -J= ^2 1 p\t k a z lfk + h.c.. 


Here ^ = (fj, fi) and n/’k = (4 t ,c_ k +) are Nambu 
spinors for the impurity and the bath, a X}Z are Pauli 
matrices in the Nambu space, e/ is a measure of the 
deviation from particle-hole symmetry on the impurity, 
U is the Hubbard repulsion, Ak is the pairing func¬ 
tion of momentum k in the bath, and finally tk is the 
momentum-dependent coupling amplitude between the 
impurity and the bath. We assume Ak = A and fk = t 
as in Ref. [12, 13]. However, we assume a more general 
normal state density of states (DOS) p(w) = po(l + kw) 
for u € [—Do, Do] and po = 1/2 Dq. Notice that k is a 
measure of the asymmetry in the DOS. Henceforth we 
take Do = 1 as the unit of energy. Integrating out the 



FIG. 2. The first excitation energy £ as a function of U for 
A = 0.01 (gray solid line) and A = 0.1 (black solid line). 
The dashed line is the result for A = 0.01 but with proper 
accounting of the under-estimation in Ref. [13] for comparison. 
See the text for details. 


superconducting bath, we get £(z) and subsequently 

= r [M g o + kw £< j z - Ac^signM] ^2) 

7Tf 

where T = np^t 2 , <to is the 2x2 identity matrix, and 
we use £ = \foj 2 — A 2 for brevity. Henceforth we use 
W(u) = 6(u)8( 1 — u) as a window function. 

First we consider k = 0 and e/ = 0, a case studied in 
Ref. [12, 13], but we proceed using the present scheme for 
comparison. We map the system to a Wilson chain as 
described above, and we perform NRG iterations there¬ 
after.We monitor the energy E of the lowest excited state 
on top of the many-body ground state as a function of U. 
The results are shown in Fig. 2 for A = 0.01 (gray solid 
line) and A = 0.1 (black solid line). The NRG results 
in the free limit U = 0 is in excellent agreement with 
the analytical results (arrows) from the Green’s function 
method. With increasing U, we see, in either case of A, 
a cusp at which E = 0. It is known that this corresponds 
to a singlet-doublet transition of the many-body ground 
state. Our transition point for A = 0.01 (gray solid line) 
is slightly larger than that in Ref. [13]. This is however 
not an inconsistency. In fact, the discretization scheme 
in Ref. [13] under-estimates the effect of T by a factor of 
A = f in A, [3, 9] where A is a scaling factor. For 
A = 2.5 used in Ref.[13], A ~ 1.069. Thus for a fair 
comparison, we use T' = T/A in our calculation, and 
present the data at T instead. The result is shown as the 
dashed line in Fig. 2, which is now in agreement with that 
in Ref. [13]. 

Next we consider both n = 0 and n = 1, and ask 
how k influences the ground state. The case of n ^ 0 
is beyond previous reach, but no excess difficulty arises 
in our strategy. For our purpose we calculate the static 
magnetic susceptibility Xi mp (T) = Xtot (T) - Xtot( T ) at 
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FIG. 3. The phase diagram for an s-wave superconducting 
bath, (a) n = 0 and (b) K = I. In both cases 7rr = 0.5 and 
A = 0.1. The color indicates the magnitude of M 2 . At zero 
temperature, the phase boundary would be sharp and mark 
the transition from doublet to singlet ground states. 



FIG. 4. The phase diagram for a d -wave superconducting 
bath, (a) t 2 /£o = 0 and (b) £ 2/^0 = 0.1. I 11 both cases to is 
fixed by 7rr = 0.5. The color indicates the magnitude of M 2 . 


an energy scale T ~ 10 6 (the zero temperature limit). [3] 
Here, Xtot (Xtot) the total susceptibility calculated with 
(without) the impurity. An effective moment scale M 
can be defined by M 2 = Txim P {T). One expects M 2 to 
vanish if the impurity spin is screened, while M 2 = 1/4 
if a full local moment survives. In fact, they correspond 
to the singlet and doublet ground states, respectively. 
The resulting phase diagram revealed by M 2 is shown in 
Fig. 3. While the phase boundary is symmetric in ef for 
k = 0 in (a), it is boosted toward the right for At = 1 in 
(b). This can be understood as follows. Given the above 
2?(w), it can be shown that At alone leads to a self-energy 
correction Sea z as oj —> 0, with 5e = —2acT(\/ 1 + A 2 — 
|A|)/7T. This means the effective impurity level is shifted 
roughly by 5e. Therefore, if T/U = g(ef/U) describes 
the phase boundary for At = 0, it becomes T /U = g[(ej + 
de)/U] for At ^ 0. This explains qualitatively the behavior 
of the phase boundary in Fig. 3(b). 

Quantum impurity in a d-wave superconducting bath: 
We now consider a d-wave superconducting bath. For 
simplicity, we use the same normal band with At = 0 as 
before. In case of constant fk, 2?(w) is actually diago¬ 
nal, the effect of which is equivalent to that of a nodal 
normal metal that can not screen the impurity spin at 
Cf = 0. [13] In this respect, it is interesting to ask how an 
off-diagonal (channel-mixing) term in T>(uj) would occur 
and how it would affect the fate of impurity spin at zero 
temperature. This is an issue not yet addressed. 

In fact, to induce channel-mixing in our case, all that 
we need is a fk asymmetric under the point group. For 
example, in a lattice model of the bath, if the impurity is 
coupled, via hopping integral f, only to two sites at r = 
ztx/2 on a nearest-neighbor bond, we would have fk = 
2tcos(k x /2) (in a suitable gauge). This can be resolved 
into A\ g and B\ g lattice harmonics. Thus in general, 
we may assume fk —> t# = ^Zitie 1 in the continuum 
limit, with </> the azimuthal angle of k, and f_j = £/ 
required by time-reversal symmetry. In the same limit 
we write Ak —> = Acos(2</>). For concreteness, we 


assume f^ = to + 2t2 cos(2<j>), and we fix A = 0.1 in the 
calculations. 

The d-wave superconducting bath leads to 
dc/>r|4| 2 (u;ero - A ^ cj x ) sign (w) 

27 t WoH, 

where T = 7rp 0 to and £<t> = \J U>2 ~~ A 2 . We make the res¬ 
olution T>(ui) = do(u))cro + d x (ui)a x to find, for |u;| <C 1, 
do(w) ~ (tp + 2t|)|w|/A and d x (u) ~ tot 2 w|w|/A 2 . Thus 
t 2 leads to a channel-mixing component d x (ui). We per¬ 
form NRG using the present strategy. Fig. 4 shows the 
resulting phase diagram revealed by M 2 , with 7rr = 0.5 
fixed. There is a marked difference when e/ = 0, the 
particle-hole symmetric point for the impurity: while a 
full local moment persists for f 2 /to = 0 in (a), we do find 
a transition from the local moment phase to complete 
screening as T/U increases for f 2 /fo = 0.1 in (b), simi¬ 
larly to the case of s-wave bath in Fig. 3(a). Indeed, the 
behavior of d x (u ;) is such that it vanishes for small |w| 
but is finite above the gap edge, resembling qualitatively 
that from an s-wave bath. This indicates the essential 
role played by an asymmetric fk- 

Summary and remarks: We developed a versatile 
scheme to map a general quantum impurity model to 
an open Wilson chain, resolving the difficulty caused by 
channel mixing. This opens the door toward broad ap¬ 
plications in protocol nonperturbative machineries, such 
as DCA and cDMFT, for strongly correlated electron 
systems. We illustrated the strategy by investigating 
the quantum phase transitions in two quantum impurity 
models with cases untouched before. 

We remark how our mapping scheme also applies to 
general Kondo impurity models. Suppose the bath is 
composed of two parts A and B , and only A is coupled 
via spin-exchange to the quantum impurity I. The idea 
is to map the effect of B on A, using our scheme, as an 
open Wilson chain C that starts from A. Then I ® C 
forms the suitable system for NRG iterations. 
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